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Abstract 

We examine the implications of singlet-doublet Higgs mixing on the properties of a Standard 
Model (SM)-like Higgs boson within the Peccei-Quinn invariant extension of the NMSSM (PQ- 
NMSSM). The SM singlet added to the Higgs sector connects the PQ and visible sectors through 
a PQ-invariant non-renormalizable Kahler potential term, making the model free from the tadpole 
and domain-wall problems. For the case that the lightest Higgs boson is dominated by the singlet 
scalar, the Higgs mixing increases the mass of a SM-like Higgs boson while reducing its signal 
rate at collider experiments compared to the SM case. The Higgs mixing is important also in 
the region of parameter space where the NMSSM contribution to the Higgs mass is small, but its 
size is limited by the experimental constraints on the singlet-like Higgs boson and on the lightest 
neutralino constituted mainly by the singlino whose Majorana mass term is forbidden by the PQ 
symmetry. Nonetheless the Higgs mixing can increase the SM-like Higgs boson mass by a few GeV 
or more even when the Higgs signal rate is close to the SM prediction, and thus may be crucial for 
achieving a 125 GeV Higgs mass, as hinted by the recent ATLAS and CMS data. Such an effect 
can reduce the role of stop mixing. 
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I. INTRODUCTION 



A natural way to explain the smallness of the higgsino mass parameter in the minimal 
supersymmetric standard model (MSSM) is to promote it to a SM singlet S so that the 
superpotential includes a coupling XSH u H d . This leads to the next-to-MSSM (NMSSM) l|, 
which however generally suffers from the tadpole (2, 3] and domain- wall {4, 5| problems once 
one adds self-interactions of S to avoid a visible axion. In particular, the tadpole problem 
makes it difficult to embed the NMSSM into a grand unified theory. 

We have recently pointed out in Ref. [6j that all the problems arising due to the SM 
singlet S can be avoided when the NMSSM is extended to incorporate the Peccei-Quinn 
(PQ) symmetry solving the strong CP problem [7, 8|. The PQ symmetry protects S from 
acquiring large tadpoles. Furthermore, the domain-wall problem is resolved by introducing 
an appropriate number of PQ messengers or considering the situation where the saxion is 
displaced far from the origin during inflation. In the PQ-invariant extension of the NMSSM 
(PQ-NMSSM), S plays the role of a messenger that connects the PQ sector with the visible 
sector. This is done through its non-renormalizable coupling with a SM singlet responsible 
for spontaneous PQ symmetry breaking at a scale much higher than the electroweak scale. 
Such non-renormalizable coupling generates a small effective tadpole for S. As a result, the 
electroweak scale originates from the SUSY breaking scale and the axion decay constant. 
Another interesting property is the existence of a relatively light neutralino dominated by 
the singlino. 

The inclusion of S modifies the Higgs and neutralino sectors, and opens the possibility 
to have singlet-doublet Higgs mixing. Such mixing increases the mass of the SM-like Higgs 
boson when the lightest Higgs boson is singlet-like 9|. It should be noted that the mixing 
effect is important even in the region with large tan/3 and small A, where the NMSSM 
tree-level contribution to the Higgs quartic coupling and thus to the Higgs boson mass is 
negligible. The PQ-NMSSM, in which quadratic and cubic terms in S are forbidden by the 
PQ symmetry, includes a light singlino-like neutralino. This results in that the Higgs sector 
has a different phenomenology from other NMSSM models because, if kinematically allowed, 
the Higgs bosons and the Z boson invisibly decay into a pair of neutralinos through couplings 
proportional to A. Furthermore, the Higgs quartic coupling can receive additional sizable 
radiative corrections involving the Yukawa interaction, Higgs-higgsino-singlino [6j. In this 
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paper, we examine the implications of the Higgs mixing on the properties of a SM-like Higgs 
boson by combining the experimental constraints placed on the singlino-like neutralino and 
the lightest Higgs boson. 

The recent ATLAS and CMS data hint the existence of a SM-like Higgs boson with 
mass around 125 GeV 10j. To account for this, one needs a particular mechanism of SUSY 



breaking giving large stop mixing or an extension of the MSSM for superparticles having 



masses around a TeV as suggested by the gauge hierarchy problem [Uj. 1 Thus, the singlet- 
doublet Higgs mixing may be crucial for achieving a 125 GeV Higgs mass without invoking 
large stop mixing. The Higgs mixing needs to be small in order for the Higgs signal rate at 
collider experiments to be close to the SM prediction. Nonetheless, we find that it can still 
increase the Higgs boson mass by more than a few GeV while avoiding the experimental 
constraints. 

In the next section, we briefly discuss the properties of the PQ-NMSSM and the effects 
of the singlet-doublet Higgs mixing. To see the impacts of the Higgs mixing, we construct 
a low energy effective theory assuming that the MSSM superparticles except the higgsinos 
have masses around a TeV. In section 3, we explore the properties of the singlino-like light 
neutralino, on which the LEP experiments put constraints. Finally, in section 4, we discuss 
how large the Higgs mixing can be in the PQ-NMSSM and how much it can increase the 
mass of a SM-like Higgs boson. Section 5 is devoted to discussion and conclusions. 



II. PQ-INVARIANT EXTENSION OF THE NMSSM 

The PQ-NMSSM is obtained by extending the NMSSM to incorporate the PQ solution 
to the strong CP problem, and the SM singlet S added to the Higgs sector connects the PQ 
and visible sectors through the operators, 

d46K Mp~i S + J d2QXSHuHd > W 

where the PQ symmetry is assumed to be spontaneously broken by another SM singlet X at 
a scale much higher than the electroweak scale. The k term induces an effective tadpole for 
S at ~ F^/Mpi after PQ symmetry breaking and drives S to acquire a vacuum expectation 



See Refs. 
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17[ for a recent discussion of singlet extensions of the MSSM. 



3 



value. Here F a is the axion decay constant, and we drop a PQ-invariant superpotential term 
X 2 H u Hd/Mpi since it can always be absorbed into the k term by redefining S. It should 
be noted that, protected by the PQ symmetry, S does not have a large tadpole. 2 On the 
other hand, as in other NMSSM models, the A term replaces a supersymmetric \i term of 
the MSSM. Having non- negligible coupling A to the Higgs doublets, S modifies the Higgs 
and neutralino sectors in a significant way. 



A. Low energy effective action 

The scalar potential for the extended Higgs sector depends on the effective tadpole term 
of S and the SUSY breaking parameters of S and H U}d , 

V = g -^^(\H u f-\H d \r + ^\HiH d \*+\\H u H d ^ 
o Z 

+ m% u \H u \ 2 + m 2 Hd \H d \ 2 + m 2 s \S\ 2 + (A x XSH u H d - B K m 2 S + h.c.) , (2) 



where we take the same notation as in our previous paper [6(, with the exception that we 
will use B and \i instead of B cS and /i e ff, respectively. It should be noted that m 2 , A\ and 
B K are set by the SUSY breaking scale m so f t , while Wq ~ n m so { t (\X \ ) 2 / Mpi is determined by 
the axion decay constant. Assuming for simplicity that the phases of A\ and B K are aligned, 
one can always make A\, B K , A and real positive by rotating the phases of S and H u d . 
We take such a field basis throughout this paper. Using the extremum conditions, one 
can replace (B K ,m 2 Hu ,rn 2 H by (B,fM, tan/3, v), where tan/3 is the ratio between Higgs 
vacuum expectation values with v ~ 174 GeV, and B is the soft parameter associated with 
the effective \i term, \i = X(\S\). Then, the Higgs sector is parameterized by 

(A,m|, A x , B, /I, tan P), (3) 

in terms of which m 2 , is written m 2 , = (B — A\)fi/\, implying that B > A\ in the model. 
See also the appendix for the relations between other SUSY breaking terms and the above 
parameters. 



The tadpole problem can be avoided also by imposing a discrete R symmetry as discussed in Refs. 



18- 



21j . Then, the SM singlet S added to the Higgs sector only has a superpotential term SH u Hd and a 
small tadpole term induced after SUSY breaking. See also Ref. [22j, where the NMSSM is extended to 
incorporate the PQ symmetry and local R symmetry. 
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The region of parameter space that we will investigate is such that the singlet scalar 
is relatively light, around the electroweak scale. To see the impacts of the singlet-doublet 
Higgs mixing, we further assume that heavy Higgs scalars and MSSM superparticles except 
the higgsinos obtain TeV masses. Integrating out the heavy fields, one obtains a low energy 
effective theory, which reduces to the SM with additional fields: the singlet complex scalar 
S, singlino S and higgsinos H Ujd . The effective scalar interactions read 

V eS = ?f(\H\ 2 ~ v 2 ) 2 +(^ + v 2 ) \XS -rf + {\H\ 2 - v 2 )\\S - 



+ - (2/i - A x sin 2/3)|(|^T| 2 - v 2 )(\S - //) + h.c.} 



sin 2/3 



/^ x sin 2 /3-/ mix cos 2 /3|V| 2 , (4) 



2Bn 

where H denotes the light Higgs doublet scalar in the decoupling limit. The last term in 
the potential results from the tree-level exchange of the heavy Higgs scalars having mass 
~ (2 J B / u/sin2/3) 1 / 2 > m w : 

/ mix = A\(\S + S ^(g 2 + - 2A 2 )(|tf | 2 - v% (5) 

which vanishes at the vacuum. The derivation of / m ; x is presented in the appendix. As in the 
conventional NMSSM, the Higgs quartic coupling is given by A#(m so ft) = 9 + 4 9 cos 2 2/3 + 
Y sin 2 2(3 at the tree-level, and receives threshold corrections coming from stop loops. The 
Yukawa interactions relevant to our discussion are 

- C cS = y t i R QH c + y b b R QH + y' s SH u H d + y' u HH u S + y' d H c H d S + h.c, (6) 

where H c = w^H* ■, and the couplings are given by y' s = A, y' u = Acos/3 and y' d = Asin/3 at 
the scale m so ft. 



B. Higgs mixing 

In this subsection, we briefly review the effects of Higgs mixing in an extension of the SM 
where the SM Higgs boson mixes with a SM singlet scalar. The low energy effective theory 
(jlj) belongs to this class of models. The CP-even neutral scalars h and s, coming from H 
and S respectively, compose the mass eigenstates Hip- 

H\ = s cos 9 — h sin 6, 

H 2 = ssin9 + hcos9, (7) 
5 



which have mass m 2 Hi = M 2 S — M\ s tan 9 and m^ 2 = M 2 h + M\ s tan 9 for the Higgs mixing 
angle fixed by 

2M 2 

where M 2 - = (didjV) is the mass squared matrix element for (h,s), and should satisfy the 
stability conditions M 2 hjSS > and (Ml) 2 < M 2 h M 2 s . 

Let us consider the case where the lightest Higgs boson Hi originates mainly from the 
singlet s. Then, the mass of a SM-like Higgs H 2 



m 



2 H2 = M 2 hh + (M 2 hh -m 2 Hi )t^ 2 9 (9) 



receives a positive contribution from the mixing. To see how much it can increase m# 2 , one 
should take into account the experimental results on the Higgs search. For a singlet-like Hi 
lighter than 114 GeV, LEP places stringent constraints on e + e~ — > ZHi — > Zbb, i.e. on the 



effective coupling [23 1 



R% = Bv(Hi bb) = Bv(Hi bb) sin 2 9, (10) 

9zZh 

where Bi(Hi — > bb) is the branching ratio of the corresponding process. On the other hand, 
the production of H2 at hadron colliders proceeds essentially through the same processes as 
those in the SM. However, the Higgs mixing modifies the discovery reach: 



R s ™ = -^r Br(# 2 -> SM) = Bi(H 2 -)• SM) cos 2 9, (11) 



which provides the signal rate for the decays of H2 to SM particles in comparison with the 
SM case. Here <j(H 2 ) is the cross section for the H 2 production, while crg M (/i) is the Higgs 
production cross section in the SM. If the Higgs signal rate is measured to be close to the 
SM prediction, only small mixing is allowed. 

The contribution to m# 2 from the singlet-doublet Higgs mixing is estimated by 

Am H2 =m H2 -M hh = — tan B + O^—), (12) 

in which tan 9 should be small enough to satisfy the LEP limits on , and will be further 
constrained if one requires the H 2 signal rate to be close to what the SM predicts. For H 2 
having mass around 125 GeV, if for instance Rj^ > 0.7 is imposed, the contribution to 
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mn 2 from the Higgs mixing can be as large as about 7 GeV when H\ has mass in the range 
between about 90 and 100 GeV. Here we have used that the LEP limits on require 
sin 9 less than about 0.25 for H% having mass in the indicated range, and that Rjjg > 0.7 
translates into sin 2 # < 0.3 assuming that H 2 dominantly decays into SM particles. Even 
if is to be more close to unity, for instance larger than 0.9, Am# 2 can be about 3 
GeV at rriHx around 90 GeV. The contribution from the Higgs mixing can thus be crucial 
for achieving a 125 GeV Higgs mass in the supersymmetric SM. In the PQ-NMSSM, the 
situation changes a bit due to invisible Higgs decays into neutralinos, but not significantly 
if one requires i?frf > 0.7. We will return to this point in section ITVl 

III. SINGLINO-LIKE NEUTRALINO LSP 

Having no Majorana mass term, the singlino obtains a small mass through mixing with 
the neutral higgsinos, and becomes the dominant component of the lightest neutralino Xv 
This opens the possibility that invisible decays of the Z and Higgs bosons into neutralinos 
are kinematically accessible. Here we assume i?-parity conservation. Then, the size of the 
higgsino component in Xi is severely constrained by the LEP limit on the invisible Z decay 
rate. Also, the Higgs production rate at collider experiments is modified. 

Let us examine the properties of the lightest neutralino before moving to discuss the 
Higgs properties. We note that the neutralino properties discussed here are independent of 
the details of the Higgs sector. In the low energy effective theory, the neutralino sector is 
comprised of two neutral higgsinos and the singlino, 



where i = (1,2,3), and we are assuming that gauginos are heavy and decouple below the 
SUSY breaking scale. The lightest neutralino couples to the SM particles through its higgsino 
component, whose size depends on e = Xv/fi and tan/3. The relevant neutralino interactions 
are 



where s and h are the CP-even neutral Higgs boson coming from S and H, respectively, and 



H° d = N 3tX l H Q u = Nuxl S = N 5iX i, 



(13) 



(14) 
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(V^) T = (Xi ) Xi) is the four- component spinor. The neutralino couplings read [11 1 



4 



A 



v/2 



{NsiNsj sin /3 + iV 4i iV 5j cos /?) + (i <-»• £ 



{NvN^-NvN^ + iH+j), 



2 cos #w 

at the scale m so f t . 

For e = Au//i <C 1, we find the neutralinos to be 



4. j ■ 



(15) 



/ X^ 

V xS/ 



V 



i _ 

1 



1+sin 2/3+cos 2/3 2 1 

4\/2 a/2 
_ 1-sin 2/3+cos 2/3 ^2 1 



1+sin 2/3-cos 2/3^2 sin/3+cos/3 \ 



4^2 



V2 



-e COS /3 



V2 



1-sin 2/3-cos 2/3 2 sin^-cos^ 
4^ 



— e sin /3 



/ho\ 

K 
\ § J 



(16) 



with masses given by 



///,.<< c± e 2 yu(l - e 2 ) sin 2/3, 



m X° 2 ~ V ( 1 + 



1 - sin 2/3 



Xi 



(17) 



where we have ignored higher order terms in e and small mixing with the gauginos. Thus, 
Xi is the lightest superparticle (LSP) for small e, and becomes lighter as tan/3 grows. 

For gauginos with mass much larger than /x, the lightest chargino Xi originates mainly 
from the charged higgsino. Thus, the LEP bound on the chargino mass demands /i > 
100 GeV. In the case under consideration, the next-to-LSP (NLSP) is either x\ or the 
higgsino-like chargino Xi ■ The 1-loop correction to the chargino mass is dominated by that 
from gauge boson loops, which is about a few hundred MeV 24J. As being higgsino-like, 
the chargino dominantly decays through the interaction W~XiXi f° r A* > m w + m x°- On 
the other hand, the main decay modes of are hito Zx\ and into a Higgs boson plus x\i 
depending on the size of \x. For e C 1, the singlino-like neutralino couples very weakly to 
other superparticles, and therefore supersymmetric cascade decays will proceed first into the 
NLSP through the same interactions as in the MSSM [25]. The NLSP then promptly decays 
into the LSP. 

Now we look into the constraints placed on the neutralino sector. If x\ has a mass smaller 
than half of the Z-boson mass, the coupling g n should be sufficiently small in order not to 
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FIG. 1: Neutralino properties for tan/3 = 5 (left) and tan f3 = 10 (right). The shaded 
region is excluded by the LEP bound on the invisible Z decay rate. Red lines represent the 
contours of the cross section for neutralino pair production a(e + e~ — > X2X1) a t \/s = 209 
GeV in the (/i, A) plane. In the right side of the dashed red line, m x o + m^o is larger than 
209 GeV. We also show the contours for the mass of Xi m t ne dashed blue line. 



exceed the experimental bound on the invisible decay rate of Z [261 . 127j . This requires 



96/t 



m z 



4m 2 o 
xi 

2 

VTir, 



3/2 



< 2 MeV, 



which translates into 



0.3 



cos 2 2/3 <1.5, 



(18) 



(19) 



for e «C 1. The above constraint becomes strong at large tan/3, where we thus need a small 
A or large /i. For fi > 100 GeV as required by the chargino mass bound, the decay Z — > X2X1 
is kinematically forbidden. 



LEP has placed constraints also on the neutralino production rate 28fl. Having higgsino 
components, neutralinos are produced at e + e~ colliders through Z exchange in the s channel. 
For fi much smaller than the neutral gaugino masses, production processes through the 
exchange of selectron in the t and u channels are negligible because of a small coupling for 
the higgsino-electron-selection Yukawa interaction. These processes are further suppressed 
when the sleptons are heavy. The LEP bound on e + e~ — > X2X1 requires 

^(e + e- x MX°2 Z X ° X ) < 70 fb, (20) 
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putting constraint on the size of #11,12 • If A* is large enough to kinematically allow the decays 
of X2 i n t° a Higgs boson and x\i the branching fraction for X2 ~~ %Xi I s reduced. 

Fig. [T] illustrates the properties of the neutralino sector, which depend on (A,/i,tan/3). 
The shaded region, where r z _ s>x o x o is larger than 2 MeV, is excluded by the LEP bound on the 
invisible Z decay rate. In the allowed region with tan /3 > 5 and 100 GeV < < 200 GeV, 
we find the following properties. The lightest neutralino has mass less than a few GeV, and 
the neutralino production cross section a(e + e~ — > X2X1) I s l ess than 70 fb in the non-shaded 
region. The constraint (120]) is thus evaded in that region regardless of how strong the process 
X2 into a Higgs boson and Xi is- On the other hand, if one takes tan (3 < 3 and A > 0.6, it 
is possible to render Xi to have a mass larger than 

We close this section by briefly discussing the LSP relic abundance in the case with 
e C 1 and fi around a few hundred GeV. 3 Since Xi is then much lighter than the Z and 
Higgs bosons and interacts very weakly with them, its annihilation processes mediated by 
the exchange of Z or Higgs bosons are too weak to avoid overdosing the universe. In order 
to resolve this cosmological difficulty, one may consider a superparticle lighter than Xi, such 
as the axino or the gravitino, into which Xi decays. For those two, the relevant interactions 
include 

C = Caaa^xi^a) + C 3/2 Gy^f^a) + h.c, (21) 
with couplings given by 

n ( K xs) « 



C3/2 ~ ZT^TT- > ( 22 ) 



m v o 1 
a-i 

m 3/2 Mpi 

where a is the axino, and a is the axion. The C a interaction comes from the k term, while 
the other is the effective interaction for the goldstino component G1/2 of the gravitino 31] 
for the case that the gravitino mass m 3 / 2 is smaller than m x o. If the axino or the gravitino is 
the LSP, Xi w iH decay with a very long lifetime. Such late decay of Xi would cause another 
cosmological problem because the produced LSP becomes a hot dark matter com pon ent, 
whose energy density is severely constrained by the CMBR and structure formation [32j, |33|. 



3 See Refs. 29|, |30( for the case with not so small e. In such case, \i obtains a mass larger than a few ten 



GeV, and its thermal relic abundance can be consistent with the observed amount of dark matter. 
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It would be possible to avoid the bound on the hot dark matter abundance for a tiny LSP 
mass. Then, the LSP, produced not by the late decay of xli an d the axion can constitute 
the cold dark matter of the universe. 

On the other hand, if i?-parity is conserved and Xi is the LSP, to avoid the LSP overpro- 
duction, one needs a sufficiently low reheating temperature or a mechanism of late entropy 
production. Indeed, in the PQ-NMSSM, it is plausible that the saxion potential energy 
dominates the universe as it has a very flat potential lifted by SUSY breaking. Let us as- 
sume such a case and examine the LSP production processes. In the absence of the k term, 
the main decay mode of the saxion a is generally a — > aa. However, the k term induces the 
interactions of cr, as well as of a, with the NMSSM particles through kinetic mixing of the 
scalar s and fermions in X and S: 

{K ^-wrA' (23) 

The saxion couplings to SM particles are ~ \m so f t /F a for ~ ra^, and are crucial for 
suppressing the branching ratio of the saxion decay into axions }34j . The LSP is produced 
by thermal and non-thermal processes. For thermally produced LSPs not to overdose the 
universe, the saxion decay temperature should be low enough. On the other hand, non- 
thermal production proceeds through the interaction (JXiXi having coupling ~ yl^Kxg) 
where y{ t ~ Ae 2 sin 2/3. Since the involved coupling is further suppressed by e 2 , it would not 
be difficult to make the produced LSP energy density small. Meanwhile, if the saxion decay 
into the axino is kinematically allowed, one should pay attention also to axinos directly 
produced by saxion decays. The axino decays via the Cg, coupling in fT2~Tl) . but mainly 
through the Higgs-a-x? 2 interactions if kinematically accessible. If lighter than the Higgs 
bosons, the axino will dominantly decay to Xi with a very long lifetime. Since the produced 
LSP then behaves like a hot dark matter, its energy density should be small. A more detailed 
analysis will be given elsewhere. 

IV. HIGGS SECTOR 

In this section, we explore quantitatively the implications of singlet-doublet Higgs mixing 
on the mass and the production rate of a SM-like Higgs boson in the case where the lightest 
Higgs boson is dominated by the singlet scalar. 
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A. Higgs properties 



Let us first examine the Higgs properties within the effective low energy theory below 
m so ft. From the Higgs potential (|IJ), it is straightforward to get 

M 2 hh = 2\ H v 2 , 

Ml = m% + A V (1 - sin 2/3 cos 2 2/3 
y 2x5 jj, 

M 2 hs = Xv(2fx- A A sin2/3). (24) 



Here radiative contributions to M? h are easily evaluated by solving the renormalization 



35 



group 
. The 



(RG) running equation of the Higgs quartic coupling X H in the effective theory 
RG running is affected by the singlet and neutralinos when they are light. On the other 
hand, the CP-odd neutral Higgs boson A originates from the singlet scalar and obtains mass, 

m 2 A = m 2 s + A V (l - ^ sin 2 $) > ( 25 ) 

which differs from M 2 S due to the contribution from / m j X . 

The main purpose of this work is to explore if the model can accommodate H2 around 
125 GeV for a singlet-like Higgs Hi having mass above about 80 GeV, below which, as we 
will see, only small mixing is allowed. In such a case, the decays of Hi 2 into SM particles are 
dominated by those into bb and VV* where V = (W, Z). Since the processes are mediated 
via the h component, one finds the decay rates to be 

Th 2 ^ sm = r S M("iH 2 ) cos 2 ^, 

IV^sm = T SM (m Hl ) sin 2 9, (26) 

in which IT sm(^ hJ is the decay rate of h — > SM obtained for a SM Higgs boson h having 
mass rriHi, in the limit of no mixing with s. In addition, because the singlino-like neutralino 
is very light, invisible Higgs decays are possible: 

k H - 2 

r Ha ^ x o x o = (y^ cos 6 + y s ij sin 6) 2 m H2 , 

r ^x°x° = -^r(^ sin ^-2/ij cos ^) 2m ^i> ( 27 ) 

for % > j, where neglecting the mass of Xi we have kf±' 2 ~ 1/2 and kf^' 2 ~ (1 — m 2 /m 2 H ) 2 
if the corresponding process is kinematically allowed. The Yukawa couplings responsible for 

12 



the Higgs decays into neutralinos read 




V2A(l + 0(e))esin2/3, y{ 



A 



(1 + C(e))e 2 sin 2/3, 



V2 



2/21 = 2(l + 0(e))(cos0-sin0), ^ 



-(l + C(e))(cos/3-sin/3)e, (28) 



at the scale m so f t . 

The neutralino Yukawa couplings to H\p show that the invisible Higgs decay H\ 2 —> XiXi 
becomes weak at large tan/3. We also note that y 21 is stronger than the bottom Yukawa 
coupling as long as A is larger than about 0.1 and tan/3 is not close to unity. Hence, if 
kinematically open, the mode Hi — » X2X11 followed by the decay of x 2 % (or a Higgs 
boson) and x\i wm dominate the Higgs decays while making Higgs searches much more 
difficult. Meanwhile, from the couplings Vi{ s , it is straightforward to see 



where ch 1 — 2 + ecot^ and ch 2 — 2 — etan# for e < 1, neglecting the masses of the final 
states. The branching ratio for the invisible Higgs decay is essentially determined by the 
above quantity, which increases as tan/3 decreases and e grows for given tan#. The invisible 
channel of Hi relaxes the LEP constraints on it, but that of H 2 reduces the Higgs signal 
rate compared to the SM prediction. 

B. Effects of Higgs mixing 

Taking into account various constraints on the model parameters, we investigate how 
much the Higgs mixing can contribute to the mass of the SM-like Higgs H 2 , which is assumed 
to be around 125 GeV as hinted by the recent ATLAS and CMS data [10]. The Higgs mixing 
should be small in order for the signal rate of H 2 decays into SM particles to be near what is 
predicted by the SM. 4 Here we focus on the region with larger than 0.7. In the analysis, 
it is convenient to use 



4 Large stop mixing can enhance the gluon-Higgs coupling mediated by squarks, altering the rate of the 
gluon fusion for the Higgs production (3^]. However, such effects are non-negligible only when stops are 
relatively light. For sfermions having masses around a TeV or higher, one can take <j{H.2) / <Jswi{h) — cos 2 
in JI]}. 




(29) 



(A, /I, tan (3, m Hl , M hh , sin 2 9) 



(30) 
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instead of the parameters (A, /i, tan/3, B, A\, nig). In the limit of A = 0, which corresponds 
to the MSSM case, large stop mixing is required to realize Mhh around 125 GeV for stops 
having a TeV mass. The NMSSM contribution to the Higgs quartic coupling improves the 
situation, but only in the low tan/3 regime if A is less than about 0.7 as required by the 
perturbativity constraint. It is thus interesting to consider the increase of m# 2 by the Higgs 
mixing. 

Before proceeding further, we summarize the region of parameter space under investiga- 
tion here: 

0.1 < A < 0.4, 120 GeV < fi < 300 GeV, 5 < tan /3 < 20, 

80GeV<m Hl < 110 GeV, 115 GeV < M hh < 120 GeV, sin 2 # < 0.3, (31) 

where the MSSM superparticles except the higgsinos are assumed to have masses around a 
TeV. It has been taken into account that the constraint (fT8|) requires a rather small A unless 
Z —7- XiXi is kinematically forbidden, and that sin 6 less than 0.3 is needed to get > 0.7. 
Furthermore, to ensure the stability of the electroweak vacuum (M^,) 2 < M^ h M^ s , one needs 
A\ ~ yutan/3 unless A is very small. This leads us to consider A\ around or above a TeV for 
at a few hundred GeV, which is consistent with the assumption of the decoupling limit 
(2Bfx/ sin 2/3) x / 2 m w because B > A\ in the model. 

Let us explain a bit more on the above parameter range. We consider [i larger than about 
120 GeV so that H 2 — > X2X1 is kinematically closed, which would otherwise dominate the 
Higgs decay and reduce the identification capability of H 2 signals at collider experiments. 
For fi > 120 GeV and A < 0.4, Xi is lighter than about 10 GeV, and the invisible channel 
H 2 becomes strong at low tan /3 while making Higgs searches difficult. To avoid such 

a situation, we consider tan /3 larger than 5. Finally, we take tan# > 0, which is the case 
with 2/i > A\ sin 2/3, because in order to make the Higgs mixing effect sizable it is preferred 
that the branching fraction of the invisible decay is large for Hi but small for H 2 as much 
as possible. This requires tan# > as can be seen from (I29|) . 

In Fig. [21 we illustrate the effects of the singlet- doublet Higgs mixing for two cases. For 
the case with (A, /i, tan/3, M hh ) = (0.2, 180 GeV, 10, 120 GeV), the neutralinos have masses 
m x o ~ 182.7 GeV and m x o ~ 1.3 GeV, and the invisible Z decay width into neutralinos is 
r , ^->x x? ~ 0.21 MeV. In this case, the Higgs mixing can increase the mass of H 2 upto about 
7 GeV in the range R s ^ > 0.7 for H x having mass around 95 GeV. On the other hand, 
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FIG. 2: Higgs mixing in two cases: (A, /J,, tan /3, Mhh) = (0.2, 180 GeV, 10, 120 GeV) for 
the left panel, while (A, fx, tan /3, M hh ) = (0.1, 130 GeV, 15, 122 GeV) for the right one. We 
plot the constant contours of Am// 2 = mn 2 — Mhh in the blue line, and also those of 
the mass of the CP-odd Higgs boson in the dotted greed line on the (mn 1 , sin 2 9) plane. 
The dashed black line represents the H% signal rate R^, compared to the SM case. The 
shaded region is excluded by the LEP limits on the coupling of Hi to 66. This constraint 
is a bit relaxed due to the invisible mode H\ — > XiXi- For the comparison, we show the 
bound on R*jj in the dashed red line, assuming that Hi decays only into 66. 



for the other case with (A, /i, tan/3, Mhh) = (0.1, 130 GeV, 15, 122 GeV), we have m x o ~ 131 
GeV and m x o ~ 0.3 GeV. The invisible Z decay width into neutralinos is r^^o ~ 0.05 
MeV. For this case, Am# 2 can be as large as 3 GeV at rriH 1 around 90 GeV even when one 
requires larger than 0.9. Therefore the Higgs mixing can considerably reduce the role 
of stop mixing in achieving a 125 GeV Higgs boson mass. Meanwhile, in both cases, the 
branching ratio for H 2 — > X1X1 is smaller than 0.1, and thus below the detectable level at 



the LHC 
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39). 



In the low energy effective theory below m so f t , M h h is determined by the relation M\ h = 
2\hv 2 for \h renormalized at the electroweak scale. The RG running of A# is affected by 
the higgsinos for /1 less than m so ft, but only slightly for small A. Note also that M^ s and 
M^ s explicitly depend on the SUSY breaking parameters. The left panel of Fig. [3] shows 
the value of Mhh for different values of the stop mixing parameter X t = (A t — fi cot (3)/m so f t 
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FIG. 3: Higgs boson mass for the case with m so ft = 1.5 TeV, [i = 180 GeV and A = 0.2. 
The left panel shows Mhh for Xt = 0, 1, 2, v6 from the below, respectively. The case with 
Xt = y/E corresponds to the maximal stop mixing. In the right panel, the blue band is 
m,H 2 obtained for Xt = 1, mui = 95 GeV and 0.18 < sin 2 8 < 0.26, for which i?f^ lies in 
the range between 0.73 and 0.8. 



in the case with m so f t = 1.5 TeV, \x = 180 GeV and A = 0.2. 5 Here A t is the trilinear soft 
parameter for H u IrQl-, and the contribution from stop loops is maximized at X t = a/6. As 
one can see, Mhh around 125 GeV requires the almost maximal stop mixing. The situation 
does not change much at tan/3 > 5 if A is less than about 0.7 as required for the interaction 
SH u Hd to remain perturbative upto Mgut ~ 10 16 GeV. In the right panel, the blue band 
represents the Higgs boson mass m H2 for the case with X t = 1 and the same values of m so f t , 
\i and A. Here we have taken = 95 GeV and 0.18 < sin 2 6 < 0.26, which are obtained 
for ms ~ 90 GeV and Ax ~ /itan ~ (a TeV). In the indicated range of model parameters, 
has a value between 0.73 and 0.8. If one requires more close to unity, the amount 
of mass increase by the mixing will be lowered. 



5 For simplicity, we use (pMf to calculate Ml h with A# evaluated by the use of one-loop RG equations, 
where we assume the squarks/sleptons and gauginos to have a universal mass m so ft- This is sufficient for 
the purpose of our discussion on how much the singlet-doublet Higgs mixing can increase the mass of the 
SM-like Higgs boson. Note that higher- loop corrections to Mhh, which involve the strong gauge coupling 
and the Yukawa couplings of the third generation fermions, can induce a shift of a few GeV on the Higgs 
boson mass 40 1. 
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On the other hand, the mass of the CP-odd Higgs boson A reads 

< = Ml - AvJL sin 3 2/3 = m\ + m\ 2 - M* hh - AvJL sin 3 2 ^ (32) 



where we have used the relations (124j) and f[25|) . In the parameter space under consideration, 
if the singlet- doublet mixing yields Am# 2 of a few GeV or larger, A has mass, 

m A ~m Hi + — ^ Am H2 - \0(1)\ x — T -, (33) 

m Hl m Hl tarn p 

and thus is heavier than Hi by about Am# 2 , but always lighter than H 2 for < Mhh- 
Fig. [2] shows the dependence of on and the singlet-doublet mixing. Note that H 2 
with mass around 125 GeV is forbidden to decay on-shell into H\H\ or AA for 80 GeV < 
rriH 1 < 110 GeV, where Amj? 2 can be sizable. Meanwhile, in the decoupling regime, A is 
mostly singlet-like and contains only a small doublet component fixed by the mixing angle, 

(34) 

to which all the couplings of A to the SM fermions and gauge bosons are proportional. 
The Abb coupling ~ (y&tan/3)0' is at most about Xyb for B > A\, but can be similar to 
or larger than the AxiXi coupling ~ Ae 2 sin 2/3 in the parameter space fl3TT) . There are 
no AZZ and AWW couplings at tree-level, while the ZHiA (i = 1,2) coupling is highly 
suppressed by small factors since it is generated from the Zh- (doublet CP-odd Higgs boson) 
interaction whose coupling itself vanishes in the decoupling limit. One thus finds that A 
mainly decays into bb and X1X1, but evades the LEP constraints from e + e~ —> ZA — >■ Zbb 
due to the suppressed production cross section. The process e + e~ — > Z* — >• H\A (and H 2 A 
if kinematically allowed) also produces A, but only with small rates due to the suppressed 
ZH\A coupling. In addition, the decay of the SM-like Higgs boson via Hi — > A A* — > 46 is 
suppressed compared to H 2 —> ZZ* — > 46 because the Abb coupling is smaller than Xy b and 
if 2 couples to AA with coupling ~ X 2 v . 

Finally, we discuss the effect of mixing between the doublet Higgs bosons, which we have 
neglected since it is small when the heavier doublet Higgs obtains a mass much larger than 
the electroweak scale. In the MSSM, the light CP-even Higgs boson is composed out of H® d , 
h = (-Re(H^) sin a + Re(if°) cos«)/a/2, and thus couples to the SM particles with 

g ht i _ cos a g hb ~ b sing g hZZ _ g hWW _ . 

g% sin/3' g% cos/3' g^ z g™ w ^ [6b) 
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where gf M denotes the Higgs coupling in the SM case. The production of a SM-like Higgs 
boson at the LHC is dominated by the gluon-gluon fusion, to which the top quark loop gives 
the dominant contribution. Using this property, one can estimate the signal rate for each 
decay channel of H 2 = s sin 9 + h cos 9 compared to the SM prediction: 



a{H 2 ) Br(H 2 ->• i) Br (# 2 -> i) /cos a 



W = -' L " ' " 1 ' " ~ ' L '" " ' ' ( ) cos 2 6» (3Q) 

(h)Bi(h ->• z)| SM Br(/i ->■ i)| SM V sin/3 y 

for small values of 9 and 5 = (a — f3 + n/2). Here the quantities with the subscript, SM, refer 
to those for the SM case. For a Higgs boson at 125 GeV, the SM predicts Br(/i — > bb) ~ 0.58, 
Br(h WW*) ~ 0.22, Br(/i -»■ ~ 0.09 and Br(/i -»■ 77) ~ 2.3 x 10~ 3 with the total 



decay width I\ ~ 4.03 MeV and the production cross section a(pp — > h) ~ 15.3 pb 41]. 
Hence, for m# 2 = 125 GeV and a moderately large tan/3, one finds 

R% « (1 -0.6 5tan/3) (l -Br(# 2 ^x?X?))cos 2 #, 
i?^* « (l + 1.4 5tan/3) (l-Br(i7 2 ^x?X?))cos 2 0, 

i?^ 7 2 « (1 + 1.45 tan 0) (l - Br(# 2 -> cos 2 0, (37) 

assuming that the Higgs invisible decay is weak and the mixing angles 9 and 5 are small. 
Though being a naive estimation, the above shows the general property of the NMSSM that, 
if 5 is positive, the mixing between doublet Higgs bosons enhances the signal rate for 77 and 



WW* while reducing that for bb 42J. In the case that the lightest Higgs boson is singlet-like, 



5 is determined mainly by the mass mixing between doublet Higgs bosons h and h', 



- ml -AV . nn . in 2 (mi- AV) 1 

<5~^n : — sin 2/3 sin 4/3 sa , '- =— , 38 

4 5/i 5/i tan 2 /3' 



for m 2 , ~ 1\B{i\l sin 2/3 3> m, 2 ^, because the contribution from the s 2 -term to 5 requires both 
s-h and s-h' mixings and is further suppressed by the mass ratio, m 2 s /{2\B^,\/ sin 2/3). The 
above indicates that the NMSSM contribution can make 5 positive when X 2 v 2 > m 2 z , i.e. 
when A is larger than about 0.52. Such effect can be sizable if h! has a mass not far above 
the weak scale. Meanwhile, doublet mixing decreases the mass of H 2 by a small amount, 

(ml - X 2 v 2 ) sin 4/3 . . . 

Am H2 s ^ ~-^- A " 5, 39 

for \8\ <C 1. However, the tree-level relation M^Jtree = m 2 z + (\ 2 v 2 — m|) sin 2 2/3 implies that 
Mhh is larger than mz at the tree- level for X 2 v 2 > m z . Thus, when 5 is positive, the decrease 
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of the Higgs mass by the doublet mixing can be compensated by the NMSSM contribution 
to the Higgs quartic coupling. 

Let us now see the situation in the PQ-NMSSM, where A > 0.52 requires /x > 260 GeV 
for a moderately large tan /3 in order to satisfy the constraint on the invisible Z decay 
rate. Provided arg(A A ) = arg(5 K ) or arg(A A ) = arg(B K ) ± ir, all the couplings in the 
Higgs potential (T5]) can be made real without loss of generality through an appropriate field 
redefinition. In the discussion so far, we have for simplicity assumed axg(A\) = axg(B K ), for 
which CP is not spontaneously broken in the Higgs sector because one can always take a 
basis where all the couplings involved are real and positive. Also note that, in this case, B 
is larger than A\ as follows from the minimization condition, and A\ ~ /x tan /3 is needed to 
avoid too large singlet-doublet mixing. Hence, when Hi is singlet-like, the heavier doublet 
Higgs will be very heavy for A > 0.52 and /x > 260 GeV, making the doublet mixing effect 
small. However, the situation changes for &rg(A\) = axg(B K ) ±7r. In this case, depending on 
the model parameters, the minimum of the Higgs potential can still lie on a point preserving 
CP. Then, since B is not necessarily larger than A\, it is possible for the heavier doublet 
Higgs to get a relatively small mass so that the doublet mixing can yield a sizable positive 
8. 

V. CONCLUSIONS 

We have examined the implications of singlet-doublet Higgs mixing on the properties of 
a SM-like Higgs boson in the PQ-NMSSM, which incorporates the PQ symmetry solving 
the strong CP problem and does not suffer from the tadpole and domain-wall problems. For 
the case where the lightest Higgs boson is dominated by the singlet scalar, the Higgs mixing 
increases the mass of the SM-like Higgs boson while reducing its couplings to SM particles. 
Such mixing effect can be sizable also for large tan and small A, where the NMSSM direct 
contribution to the Higgs mass is negligible. However, the amount of mass increase is limited 
by the LEP constraints on the properties of the singlet-like Higgs boson and the lightest 
neutralino constituted mainly by the singlino. In addition, in order for the Higgs signal rate 
to be close to what the SM predicts, the mixing should be small. Combining these, we find 
that the mixing can enhance the Higgs boson mass by a few GeV or more even when the ratio 
for the Higgs signal rate compared to the SM case is to be 0.9. For the ratio around 0.7, the 
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amount of increase can be as large as about 7 GeV for the singlet-like Higgs boson around 
95 GeV. Thus, the singlet-doublet Higgs mixing may be crucial for achieving a 125 GeV 
Higgs mass, as hinted by the recent ATLAS and CMS data, within the supersymmetric SM 
with superparticles having masses around a TeV. Once the Higgs production cross section is 
measured, we will learn the bound on how much the Higgs mixing can modify the properties 
of the SM-like Higgs boson. 
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Low energy effective Higgs potential 

In this appendix, we present the relations between (A, m|, A\, B, fi, tan 8) and the SUSY 
breaking parameters involved in the Higgs potential. We also discuss how to integrate out 
the heavy Higgs scalars in the decoupling limit. 

The SUSY breaking parameters are written 

r 1 ( 2 , A tan/3 A x \ 2 2 \ 

B K = — m q +1 ttt, ) X v , 

B — A\ \ V tan 2 8 + 1 // J J 

2 B{i 2/i 2 + AV tan 2 8 - 1 / ra| , 

mHu ~ tan/3 ~ tan 2 + 1 ~ tan 2 (3 + 1 I ~2~ + ^ 



Bu 2u 2 + X 2 v 2 cot 2 8 - 1 f ml , 

m -"-J-r^TT-^fTT{-f + " )• (40) 



which are obtained by using the extremum conditions. 

Let us construct an effective Higgs potential by integrating out the heavy Higgs scalars 
in the decoupling limit (2Bfi/ sin2/3) 1 / 2 ^> mw- m the field basis H = — if^sina + H° cos a 
and H' = H d cosa + Hasina with H° = ia 2 H*, the Higgs potential reads 

V = V (S, \H\ 2 ) + f 2 \H'\ 2 - | (/a cos 2 a - f* sin 2 a)H ] H' + h.c.} + ■ ■ • , (41) 
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where fi^ are given by 

h — m 2 H cos 2 a + m 2 H sin 2 a + A 2 |S1 2 — -X(A\S + ml + h.c.) sin 2a, 

fi = ^(m 2 Hd ~m 2 H Jtan2a + X(A x S + m 2 )~^(g 2 + g /2 -2X 2 )\H\ 2 sm2a, (42) 

and the ellipsis indicates Higgs quartic terms with two or more powers of H'. It is straight- 
forward to find 

(ft) oc cos(a-/3). (43) 

Therefore, (fi) vanishes at a = (3 — tt/2, for which (H°) = v and (H') = 0. In this case, if 
2Bn/ sin 2(3 3> v 2 , the heavy Higgs doublet H' can be integrated out by solving du'V = 0: 

H' = (A* sin 2 (3-h cos 2 + ■ ■ ■ )H, (44) 

where the ellipsis in the bracket includes terms depending on S and \H | 2 , which are irrelevant 
to the low energy physics because fi vanishes at the vacuum. The relevant interaction terms 
arise only through the dependence of fi on S and \H\ 2 . Using the extremum conditions, one 
can rewrite fi, 

fi = A X (XS -ri + S ^-(g 2 + g ' 2 - 2X 2 )(\H\ 2 - v 2 ) = f mix . (45) 
Finally, substituting H' by the solution (jUJ) leads to the effective potential (jlj). 
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